Role of the spin-orbit splitting and the dynamical fluctuations in the Si(557)-Au 

surface. 
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Our ab initio calculations show that spin-orbit coupling is crucial to understand the electronic 
structure of the Si(557)-Au surface. The spin-orbit splitting produces the two one- dimensional bands 
observed in photoemission, which were previously attributed to spin-charge separation in a Luttinger 
Uquid. This spin splitting might have relevance for future device applications. We also show that the 
apparent Peierls-like transition observed in this surface by scanning tunneling microscopy is a result 
of the dynamical fluctuations of the step-edge structure, which are quenched as the temperature is 
decreased. 
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Atomic-scale "wires" created on semiconductor sur- 
faces have the potential to be the basis for nanoscale de- 
vices and to demonstrate the novel physics of one dimen- 
sion. A specially attractive route to create such struc- 
tures is provided by the use of stepped vicinal surfaces, 
which can be used as templates to deposit material creat- 
ing very regular arrays of parallel metallic chains 1]. The 
Si(557)-Au reconstruction is one of these systems. The 
terraces of this surface are perpendicular to the (111) di- 
rection and have a width of ~ 19 A. Each terrace contains 
a monatomic chain of gold atoms running parallel to the 
step-edge. For this reason the Si(557)-Au surface has 
been proposed as an experimental realization of a one- 
dimensional metal and has attracted much attention in 
recent years However, in spite of these 

efforts, the electronic structure of this system is not yet 
completely understood. 

The first angle-resolved photoemission (ARP) study 
by Segovia et a/.|^ found a spectrum dominated by a 
one-dimensional metallic band. This band was shown to 
split in two peaks near the Fermi level (Ei?), and this was 
interpreted as signature of separated charge and spin low- 
energy excitations as predicted by Tomonaga-Luttinger 
theory of the one-dimensional electron gas j^, [l^ . How- 
ever, later photoemission data seemed to discard this 
interpretation. According to Losio et al. 3] the ob- 
served splitting would correspond to two distinct proxi- 
mal bands which cross E^? at neighboring, although dif- 
ferent, positions of the surface Brillouin zone. However, 
the origin of these bands was unclear. Finally, Ahn et 
al. have recently suggested that only one of the bands 
is truly metallic and suffers a metal-insulator transition 
upon cooling. This observation was correlated with the 
temperature dependence of the scanning tunneling mi- 
croscopy (STM) images: the step-edge undergoes a peri- 
odicity doubling consistent with a Peierls-like instability. 
Thus these authors concluded that at least one of the 



two proximal bands should be associated with the atoms 
forming the step-edge. 

In contrast with the electronic structure, the geometry 
of the Si(557)-Au reconstruction seems to be quite well 
established. A detailed model was recently proposed by 
Robinson et al. 6] on the basis of X-ray diffraction data, 
which has been corroborated by first-principles density 
functional (DFT) calculations 8] . Unfortunately, the cal- 
culated band structure only presents one one-dimensional 
metallic band exhibiting a considerable dispersion and 
a width consistent with the experimental observations. 
Thus, the observed two-band photoemission spectrum re- 
mains unexplained from a theoretical point of view. 

In this Letter we present the results of DFT calcula- 
tions that satisfactorily explain the observed ARP spec- 
tra and the temperature dependence of the STM images 
of the Si(557)-Au surface. The two proximal bands ap- 
pear as a consequence of the spin-orbit (SO) splitting 
of the most dispersive surface state. The large effect 
of SO splitting on gold-derived surface states has been 
demonstrated previously by ARP experiments and cal- 
culations on Au(lll) jnlUJ. The origin and magni- 
tude of the SO splitting is similar in the present case. 
The inversion symmetry is always broken at the surface, 
thus the spin splitting of bands with no spatial degen- 
eracy becomes possible. This suggests that atomic wires 
formed by heavy atoms deposited on semiconducting sur- 
faces could be used in the fabrication of spin transis- 
tors ^3 and spin- filter devices Q . Interesting phenom- 
ena can also appear associated with the competition be- 
tween electron-electron interactions and SO coupling in 
one-dimension il5l. 11^. 

Our results also show that the temperature dependence 
of the STM images stems from the buckling of the step- 
edge, whose atoms alternate between up and down po- 
sitions. A given "up-down" configuration and the re- 
verse one are separated by a small energy barrier. At 
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FIG. 1: (a) Calculated equilibrium structure of the Si(557)- 
Au reconstruction. The corresponding electronic band struc- 
ture is shown for a calculation not including (b), and includ- 
ing (c) the spin-orbit interaction. Energies are referred to the 
Fermi level. Surface states have been marked with different 
symbols according to their main atomic character (see text). 



room temperature the step-edge fluctuates between both 
structures, and the STM only reflects the average surface 
electronic and atomic structure. Upon cooling the step- 
edge buckling can be revealed using STM, producing the 
apparent doubling of the periodicity. 

Our slabs contain four silicon bilayers, the bottom 
surface saturated with hydrogen, with a total of 115 
atoms. Due to this large supercell we use the SIESTA 
code 0, study the energy landscape and the 

structural properties of the Si(557)-Au surface. The de- 
tails of the calculations are similar to those presented in 
Refs. 1^ and for the same system. Since the present 
version of SIESTA does not allow to include the SO 
interaction, we have used the VASP code 0, to 
study its effect on the electronic structure of our relaxed 
structures. The VASP calculations utilized projected- 
augmented-wave potentials with a plane-wave cutoff of 
250 eV, which proved to be sufficient to obtain converged 
electronic band structures. 

Fig. da) shows the relaxed structure of the Si(557)-Au 
reconstruction. This structure is almost identical to the 
experimentally proposed model 0,1^. The larger, filled 
circles stand for the gold atoms occupying silicon sub- 
stitutional positions on the middle of the terraces. The 
corresponding band structure, along the direction paral- 
lel to the steps, is shown in Fig. ^b) and (c). Panel (b) 
shows the results from a non spin-polarized VASP calcu- 
lation using the local density approximation (LDA) for 
the exchange-correlation potential. This band structure 
is almost identical to that already reported in Ref. 8 using 
the SIESTA code. Several surface bands and resonances 



can be identified, all of them with negligible dispersion 
in the direction perpendicular to the steps. The different 
symbols reflect their main atomic character, as obtained 
from the projection of the charge onto non-overlapping 
spheres centered on each atom. The unoccupied band 
marked with open triangles comes from the adatoms (la- 
beled ad in Fig. ^a)), while the occupied one is related 
to the restatoms (labeled rest). In principle, every atom 
in the step-edge has a dangling-bond pointing perpen- 
dicularly to the step, which would give rise to a very 
flat half-filled band. This unstable situation leads to a 
buckling of the step-edge that doubles the unit cell and 
forms two bands marked by open squares in Fig. ^ The 
band with larger weight in the "up" (B) atoms is fully 
occupied, while the band associated with the "down" (A) 
atoms has a small occupation 21]. Notice that the step- 
edge bands have a very small dispersion. Thus they can- 
not explain the observed ARP spectra, as was recently 
suggested by Ahn et al 4] . Contrary to the initial inter- 
pretation of the ARP data 2] , none of the surface bands 
in the proximity of E^ has a clear Au 6s character. This 
is a direct consequence of the larger electron affinity of 
gold as compared to silicon: the 6s Au character appears 
several eV below Ef. The sole surface bands exhibiting 
an appreciable gold component in Fig. Q are marked with 
circles. These bands mainly come from the 3p states of 
the silicon atoms neighboring to the Au chains (atoms D 
and C). However, they also show a strong contribution 
from the 6p states of gold, so they are better assigned to 
the Si-Au bonds. We find a flat band (open circles) as- 
sociated with the Au-Sii) bonds, and a dispersive (filled 
circles) corresponding to the Au-Sic bonds. This last 
band is the only one that presents a dispersion and width 
consistent with the photoemission data. 

Fig. ^c) depicts the same band structure once the ef- 
fect of the SO interaction has been included in a non- 
collinear spin calculation. We still get a non spin- 
polarized ground state. The changes are negligible for 
most surface bands. This is expected since most of them 
are localized in regions far from the gold atoms or have 
a very small dispersion, i.e. small group velocities. How- 
ever, the dispersive Au-Sic band develops a considerable 
SO splitting. This brings the calculated band structure 
in reasonable agreement with the experimental spectra. 
The experiments 0, Q show a ~300 meV splitting near 
Ep- This splitting exhibits a linear dependence as a func- 
tion of k\\ with a ^1.2 eV A slope. This can be com- 
pared with our calculated ~200 meV splitting at E^ and 
~1.4 eV A slope. 

Even though the two SO-split bands are a robust fea- 
ture predicted for this surface, the states at the Fermi 
energy are strongly affected by the other bands shown 
in Fig. CJb) and (c). In our calculations both SO-split 
bands are metallic; however, there is a band gap just 
above Ef- This gap relates to the presence of a row of 
adatoms in the terrace, which induces an alternation of 
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FIG. 2: Total energy (a) and position of the B atom (b) as a 
function of the height of atom A. Heights are referred to that 
of atom A in the ground-state configuration. The arrows in 
panel (b) mark the position of the two energy minima. The 
energy landscape shows two wells separated by a small energy 
barrier of ~29 meV, and estimated zero-point energies of uo ^ 
5 and 9 meV. 

the Sic-Au-Sic bond angle between 101.8^ and 109.6^. 
The presence of this gap drives the surface very close to 
becoming semiconducting. The metallicity of the system 
is due to the very small partial occupation of the upper 
step-edge band that pins the position of Ei?. A descrip- 
tion of exchange and correlation beyond LDA would quite 
likely cause a widening of the step-edge gap, leading to a 
semiconducting band structure at zero temperature. 

The ground state structure and the electronic bands 
shown in Fig. ^ perfectly explain the low-temperature 
STM images of the surface U]. The step-edge exhibits 
a buckling, with A atoms lying ~0.4 A below B atoms. 
Furthermore, a high contrast between these two types 
of silicon atoms is guaranteed. Empty-state STM images 
will preferentially show A atoms, while filled-state images 
will reflect the location of B atoms. 

The origin of the temperature dependence of the STM 
images is clarified in Fig. [21 where the stability of the 
step-edge buckling respect to thermal vibrations is ana- 
lyzed using SIESTA. Fig. [2 (a) shows the energy land- 
scape as a function of the position of atom A, i.e. the 
lower atom of the step-edge in the ground-state struc- 
ture. We progressively increase its height, Za (see the 
filled symbols in Fig. [21). At each step the whole struc- 




FIG. 3: Simulated STM images at low (a) and high tempera- 
ture (b) for a +0.7 V bias voltage (empty states). The upper 
part of panel (a) corresponds to the ground state structure, 
while the lower part is obtained with the reversed step-edge 
configuration. Panel (b) combines both images. 

ture is relaxed while keeping Za as a constraint. As ex- 
pected the energy increases quadratically for small dis- 
placements from the equilibrium position {Za=0). The 
height of the B atom (Zb) remains almost unchanged 
during this process as shown in Fig. I^fb) and, although 
the magnitude of the gap decreases as Za approaches Z^, 
we can still unambiguously identify the fully and partially 
occupied step-edge bands, respectively, with the atoms 
B and A. However, for larger displacements the situation 
is very different. When Za>Zb the structure becomes 
unstable and collapses into a new configuration with an 
inverted step-edge buckling. Atom A is now situated 
~0.4 A above atom B. This sudden change is accompa- 
nied by a charge transfer from atom B to atom A, so that 
the filled step-edge band can be always assigned to the 
atom occupying the highest position. 

The reversed buckling configuration is also stable 
against small displacements (see the open symbols in 
Fig.[2j. Its energy is only ~3 meV (~5 meV) larger than 
that of the ground-state structure within the LDA (gen- 
eralized gradient approximation, GGA). For large dis- 
placements we find a transition back to the ground-state 
structure (open symbols in Fig. [21). This means that the 
roles of atom A and B can be interchanged without any 
appreciable change of the global structure of the surface 
and the total energy. This is not surprising if we consider 
that, for the structure shown in Fig. ^a), the two step- 
edge sites only become inequivalent when we consider 
their registry with the adatom row in the same terrace, 
which is located at ~12 A from the step-edge. Thus the 
energy landscape as a function of the height of the atoms 
in the step-edge is formed by two similar wells. The en- 
ergy barrier separating both structures turns out to be 
quite small, -29 meV (-25 meV) within LDA (GGA). 
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From the total energy curves and the corresponding 
displacements of all the atoms in the supercell, the zero- 
point energies are estimated to be coq ~5 meV for the 
ground-state structure, and ~9 meV for the reversed 
buckling one. This, together with the magnitude of the 
calculated barrier, shows that at low- temperatures both 
structures are indeed stable for long times. The simu- 
lated images using Tersoff-Hamann theory can be 
found in FigEl (a), the upper part of the panel corre- 
sponding to the ground-state configuration and the lower 
part to the reversed step-edge buckling. The images are 
very similar. They are dominated by two chains show- 
ing a double periodicity along the step-edge direction in 
good agreement with low-temperature images reported 
for this system 0, Q. One of the chains is the row of 
adatoms. The other corresponds to the step-edge, where 
only every other atom is visualized. 

The situation changes at higher temperatures. The 
time necessary to flip between the different step-edge 
configurations gets shorter as the temperature increases. 
Given the size of the calculated energy barrier, at room 
temperature the STM images will only show a time av- 
erage of the electronic and atomic structure of the sur- 
face. This has been modeled in Fig. |S1 (b) by averag- 
ing the STM images of the two structures. As a conse- 
quence, the double periodicity is lost and all the atoms 
in the step-edge appear with similar intensities. This 
explains the apparent periodicity doubling observed by 
Ahn et al. 4] upon cooling. Therefore, the mechanism 
behind the change of appearance of the STM images of 
the Si(557)-Au surface is similar to that proposed to ex- 
plain, for example, the 3x3^ x ^3 transition in 
Sn/Ge(lll) 23]. 

Our results also provide a plausible scenario to explain 
the observation of a metal-insulator transition in this sur- 
face 4J . The Fermi energy - and thus the occupation and 
metallicity of the SO-split bands - depends strongly upon 
the position of the upper step-edge band. Small changes 
may drive the system from metallic to insulator and con- 
versely. Furthermore, as the temperature increases the 
step-edge structure is expected to fluctuate at an increas- 
ing rate. Associated with this atomic movement there is 
an electron dynamics that can be interpreted as a succes- 
sive closing and opening of the step-edge gap or, alter- 
natively, as an increasing width of the step-edge energy 
levels. Therefore we can expect changes in the spectral 
weight at Ei? as the temperature is increased. 

In summary, we report on first-principles DFT calcula- 
tions unveiling the origin of several phenomena observed 
in the Si(557)-Au surface. The experimental electronic 
band structure is successfully reproduced, pointing out 
the important role played by the spin-orbit coupling in 
this system. From a general perspective, this raises the 
question of whether systems composed by atomic-scale 
wires of heavy atoms on semiconducting substrates can 
be used to create or transport spin-polarized currents, 



and thus be useful for future electronic devices. The tem- 
perature dependence of the STM images can be easily ex- 
plained by the presence of two structurally inequivalent 
silicon atoms in the step-edge that fluctuate at room tem- 
perature between two positions. The presence of a elec- 
tronic gap slightly above Ei? together with the step-edge 
fluctuation can also explain large changes in the bands 
near the Fermi energy as a function of temperature. 
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